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In the present work, we combine the notion of T^T-symmetry with that of super-symmetry (SUSY) 
for a prototypical case example with a complex potential that is related by SUSY to the so-called 
Poschl-Teller potential which is real. Not only are we able to identify and numerically conhrm the 
eigenvalues of the relevant problem, but we also show that the corresponding nonlinear problem, in 
the presence of an arbitrary power law nonlinearity, has an exact bright soliton solution that can 
be analytically identified and has intriguing stability properties, such as an oscillatory instability, 
which the corresponding solution of the regular nonlinear Schrodinger equation with arbitrary power 
law nonlinearity does not possess. The spectral properties and dynamical implications of this insta¬ 
bility are examined. We believe that these hndings may pave the way towards initiating a fruitful 
interplay between the notions of PT-symmetry, super-symmetric partner potentials and nonlinear 
interactions. 


I. INTRODUCTION 

In the past 15 years, there has been a tremendous growth in the number of studies of open systems bearing both 
gain and loss, motivated to a considerable degree by the study of the specially balanced T^T-symmetric dynamical 
models The original proposal of Bender and collaborators towards the study of such systems was made as 

an alternative to the postulate of Hermiticity in quantum mechanics. Yet, in the next decade, proposals aimed at 
the experimental realization of such T^T-symmetric systems found a natural “home” in the realm of optics mis]. 
Within the latter, the above theoretical proposal (due to the formal similarity of the Maxwell equations in the 
paraxial approximation and the Schrodinger equation) quickly led to a series of experiments [6]. In turn, these 
efforts motivated experiments in numerous other areas, which span, among others, the examination of T^T-symmetric 
electronic circuits lai, mechanical systems [9] and whispering-gallery microcavities pU] . 

In the same spirit, another important idea that has originally been proposed in a different setting (namely that of 
high-energy physics m) but has recently found intriguing applications in the context of wave guiding and manipula¬ 
tion in the realm of optics is that of super-symmetry (SUSY) [12]. The main idea is that from a potential with desired 
properties, one can obtain a SUSY-partner potential that will be isospectral to (i.e., possess the same spectrum as) 
the original one, with the possible exception of one eigenvalue. In fact, taking the idea one step further, starting from 
a desired ground state eigenfunction, one can design the relevant super-symmetric partner potentials in a systematic 
fashion, as discussed, e.g., in [T2|, both for continuum and even for discrete problems. In fact, more recently, the 
two ideas (of T^T-symmetry, or anyway non-hermit icity, and SUSY) have been combined to construct SUSY-partner 
complex optical potentials designed to have real spectra m- An expected application of these ideas that has also 
started to be explored (extending the spirit of corresponding studies in the T^T-symmetric setting [14]) is in using 
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SUSY transformations to achieve transparent and one-way reflectionless complex optical potentials m- 

The above works have essentially constrained the interplay of VT symmetry and SUSY at the level of linear 
potentials. Naturally, however, except for very low optical intensity, the crystals considered in the relevant applications 
bear nonlinear features, e.g., due to the Kerr effect. Hence, our focus in the present work will be to extend these 
linear ideas of VT symmetry and SUSY to a nonlinear case example. Moreover, we will select an example that blends 
two additional characteristics. On the one hand, one of our super-symmetric partners will constitute a famous and 
well-known solvable model in quantum mechanics, namely the celebrated Poschl-Teller potential [ISIIIT]. On the 
other hand, it will turn out to be the case that not only the linear but also the nonlinear variant of the problem will 
be analytically solvable, in fact for arbitrary powers of the nonlinearity, in a special limit and will naturally connect 
with the linear solutions of the potential. In what follows, in Sec. II we will first present the general theory of linear 
T^T-supersymmetric potentials. Then, in Sec. Ill we will consider the special nonlinear solutions and their asymptotic 
linear limit reduction. Numerical results will corroborate the above analytical findings and we will also explore the 
spectral and dynamical stability of the nonlinear waveforms. Finally, in section IV, we summarize our findings and 
present our conclusions. 


II. A LINEAR NON-HERMITIAN SUPERSYMMETRIC MODEL 


As is done generally in the theory of SUSY, we consider an operator A such that 

(I) 

where W is the super-potential and an operator B of the form: 

It is important to accentuate here (see also ini) that in the case of a complex super-potential kF, contrary to the 
Hermitian case of a real kF, B is not a Hermitian adjoint operator of A (hence the different symbol). Then, defining 
the potentials = kF^ ^ kF' + with — 2kF', we have that the operators 

+ ( 3 ) 

are isospectral (with the exception of the fundamental mode in the potential 1/^+^ which lacks a counterpart in 
VT), More specifically, the eigenvalues satisfy = E^^ for n > 1 (cf. also [IS]). We note in passing that the 
eigenvectors of the two cases are also related, i.e., = Au^^^ and = Bui^\ 

Now, assuming that kF = f Tig^ ^ -\-iVj ^ and that G M, we find that (cfr. [18]) 


the potentials have to satisfy the following conditions: 

vi+'> = f-g^-f + E, (4) 

= 2fg-g', (5) 

yE = f-g^ + f + E, ( 6 ) 

y-) = 2fg + g'. (7) 


The remarkable finding of the linear spectral analysis of m was that these authors, motivated by the sl(2, C) potential 
algebra were able to derive a number of special case examples of simple functional forms of complex kF’s which give 
rise to complex SUSY potentials. Perhaps the most remarkable of their examples concerns the super-potential 


W{x) 



tanh(x — c) — ibisech{x — c). 


(8) 


which gives rise (assuming hereafter without loss of generality that c = 0) to the super-symmetric partners of the 
form: 


~ sech^(x) — 2im6/sech(x) tanh(x), (9) 

f ^ j sech^(x) — 2i{m — l)6/sech(x) tanh(x). 


( 10 ) 
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FIG. 1: The blue circles represent the eigenvalues of the operator of Eq. ([^, under the T^T-symmetric potential of 

Eq. The solid lines represent the analytical predictions based on the potential’s super-symmetric partner corresponding 
to an analytically tractable Poschl-Teller potential. 


We chiefly focus hereafter on the remarkable special case of m = 1, previously considered, e.g., in [19]. The 
exceptional characteristic of this case is that it stems from a real potential which is well-known to be exactly 
solvable in the realm of elementary quantum mechanics, namely the Poschl-Teller potential [ig HI]. While its 
eigenfunctions can also be written down in an explicit form by means of hypergeometric functions, here we will 
restrict our considerations to the relevant (bound state) eigenvalues which in the context of the above example 
assume an extremely simple form as: 

Ei-^ = -\[2bi-2n-l]\ ( 11 ) 

Such bound state eigenvalues only exist when n < bj — 1/2. This, in turn, suggests that for the + superscript 
potential, it will be: i.e., all the relevant bound state eigenvalues should also emerge in the VT- 

symmetric spectrum of the potential just as they appear in the Hermitian (real) spectrum of the potential 

The only eigenvalue that will not be captured by this relation is E = —1/4; see the relevant details on the 
spectrum of 1/^+^ below. Furthermore, we expect that, when varying 6/, bound state eigenvalues will emerge as bj 
crosses 0.5, 1.5, 2.5,... in both the spectra of 

All of these conclusions are fully corroborated by the results of Fig. The spectrum of considered therein 

turns out to be real, as may be anticipated by the T^T-symmetry of the model, but more importantly, it turns out to 
be identical to that of its super-symmetric Poschl-Teller partner, as can be seen from the theoretical lines confirming 
the bifurcation of the point spectrum eigenvalues at the locations theoretically predicted. Finally, indeed, the only 
eigenvalue that is not captured is E = —1/4 which turns out to be invariant^ under variations of bj. As a final note, we 
point out that generalizations of this potential with arbitrary coefficients in both the real and the imaginary part were 
considered in m and the relevant T^T-symmetric transition threshold was identified as an inequality associating the 
real and the imaginary part prefactors. The pertinent inequality here assumes the form (6/ — 1)^ >0 and is generically 
satisfied (i.e., V 6/), as can be expected by the super-symmetric partnership of the potential with a Hermitian one 
bearing real eigenvalues for all bj. 

As a side remark, we observe that is invariant under bj O m. Further, both m and bj are arbitrary real 

numbers and not integers. Interestingly, as shown in [21], in case 6/ — m is not an integer, then remarkably, the 
eigenvalue spectrum has two branches: 

=-(m-n-1/2)2, n = 0,1,2, ...,nmax (12) 

where m — 3/2 < Umax < m — 1/2, and 

= -{hi-n-l/2f , n = 0,l,2,...,nmax (13) 

where bj -3/2 < Umax < k - 1/2. 

From this, we infer that when m = 1 and bj is not an integer, has two nodeless states (i.e. n = 0) with energy 
eigenvalues and eigenfunctions 


= -1/4, = Vsech(:c) ^ 


(14) 
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= -{hi - 1/2)2 ^ ^(2) _ g2itan-^(tanh../2) ^ 

although the latter (as per our spectral results of Fig. will only be present for bj > 1/2. Interestingly, while for 
1/2 < hi < 1, Eq = —1/4 is the ground state, but for 1 < hi < 3/2, Eq = —(6/ — 1/2)2 corresponds to the ground 
state. 


III. NONLINEAR GENERALIZATION OF THE MODEL 


We now turn to the corresponding nonlinear model which is the basis for the present analysis. Examining the case 
of the focusing nonlinearity, the operator is augmented into the nonlinear problem: 


tUt 


= - 




(16) 


The most physically relevant case is that of the cubic nonlinearity /^ = 1, corresponding to the Kerr effect, although in 
recent years, examples of higher order nonlinearities (like k = 2 and k = 3) have been experimentally realized; see, e.g., 
for a recent example [22] . The relevant nonlinear problem has been partially considered for /i: = 1 in a two-parametric 
generalization of the potential associated with [23]; see also the more recent discussions of [24l|25]. While all of 
these works were restricted to the cubic case, below we will obtain exact solutions for arbitrary nonlinearity powers. 
Moreover, we will explain, through our VT-SVSY framework the existence of nonlinear dipole (and, by extension, 
tripole etc.) solutions identified in m, emerging from the higher excited states of the underlying linear problem. It 
can be directly found that the relevant nonlinear single-hump solution for arbitrary k is of the form: 




where 


= 


4k^ 


(k + 2)2 


+ 2 ) 2 


Ah? 

[ i,? ‘'J 


(17) 


(18) 


, . , ArriKhi 

oix) = -tan 

^ kE2 


1 (tanh(|)) 


and 


(19) 




( 20 ) 


Note when m = 0 and bj = 1/2, ^ 0 and our solution reduces to the well known solution of the NLSE with 

^ 2 /^ = (1 + /^)//^2. Also note that when n = 2^ 


\ 2k, 


(6?-l)(m2-l), 


( 21 ) 


which vanishes at either 6/ = 1 or the special case m = 1. Hereafter, we again restrain consideration to the special 
case of K = 1. 


Eor bj —> bi^ci with b^ ^ , the amplitude A tends to zero and the solution (17) becomes the solution of the 

corresponding linear limit (15) by virtue of condition (20). The solution 0 exists for 5/ < bj^c when n <2 and for 
> bj^c if K, > 2. Our analytical expression only yields the trivial solution foTK = 2 as mentioned earlier. Notice 
also that when bj is fixed and k varied, the solution tends to (14) when approaching the hz = 2 limit. In addition, Eig. 
d shows the dependence of norm with respect to bj and k when the condition @ for solution @ is applied. The 
value of the norm is 


N = 


f 


\u{x,tWdx = A^B ( 1 


( 22 ) 


where it has been taken into account that A G M and B{x^y) is the Euler’s beta function. 

Horizontal “cuts” along the graph of Eig. [^ are shown in Eigs. and where the continuum tendency to the 
linear limit (dark) is shown as a variation over 6/ for /^ = 1 and /^ = 3, respectively. Apart from the analytical 
solution (1^ which collides with the nodeless solution of the linear Schrodinger equation, we have been able to find 
numerically the branch that collides with the solution with a node [n = 1 in (12)-(13)]. These solutions are the 
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FIG. 3: (Top panel) Norm of the solutions with = 1 as a function of bi. The blue line corresponds to the nodeless solution 
whereas the red line corresponds to the “dipole” branch (cf. [25]) possesing a node. The bottom left (right) panel showcases 
the modulus of the nodeless solution (solution with a node) as a function of x for different values of bj. It can clearly be seen 
that the amplitude of the solution goes to 0 as 6/ ^ 3/2 (5/ ^ 5/2). 


generalizations (for arbitrary k) of the “dipoles” of [25]. In those cases, solutions exist as long as bj < bi^c + 1 and 
their monotonicity for hi > 2 is opposite to that of the fundamental solutions analytically identified above (hence, 
the above mentioned collision). Interestingly, it is worth mentioning that the latter dipole branch is present even for 
K = 2. In the same spirit, higher order generalizations (e.g. tripoles, quadrupoles, etc.) can also be expected in the 
spirit of m, degenerating to the linear limit, respectively for bj —> bj^c + 2, 6/ bj^c + 3, etc. 

We now turn to the detailed stability analysis of the relevant soliton solutions (which was not explored systematically 
in [23] [24], but was touched upon in [25] for hi = 1). In fact, in [23] a particular case example of an evolution (cf. 
Fig. 2 therein), as well as the positivity of the Poynting vector flux led those authors to conclude that the relevant 
solutions were nonlinearly stable. However, our systematic spectral stability analysis, illustrated in Figs. and 
indicates otherwise. In particular, we use a linearization ansatz of the form 


u{x,t) = e*-®* \^uo{x) + (a(x)e^* + 6*(:c)e^*‘) 


(23) 


where uo{x) is the spatial profile of the standing wave solution of Eq. while {a(x),6(x)} and A correspond. 
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FIG. 4: (Top panel) Norm of the solutions with = 3 as a function of bj. The blue line corresponds to the nodeless solution 
whereas the red line corresponds to the solution with a node (i.e., the “dipole”). The bottom left (right) panel showcases the 
modulus of the nodeless solution (solution with a node, although the node itself disappears as bi increases) as a function of x 
for different values of bi. It can clearly be seen that the amplitude of the solution goes to 0 as 6/ ^ 5/6 {bi ^ H/^), i-e., the 
corresponding linear limit value for E = —1/9. 


respectively, to the eigenvector and eigenvalue of the linearization around the solution. The existence of eigenvalues 
with Re(A) > 0 would in this (T^T-symmetric and hence still ensuring the quartet symmetry of the relevant eigenvalues) 
context signal the presence of an instability. 

We can see in Fig. j^that indeed such an instability is present in the interval 0.56 < bj < 1.37 for the nodeless 
solutions of K = 1. Further examination of the relevant phenomenology in Fig. [^reveals the origin of the instability 
and its stark contrast with the corresponding phenomenology in the standard NFS model. In particular, the breaking 
of translational invariance (due to the presence of the potential) leads the corresponding neutral mode to exit along 
the imaginary axis of the spectral plane (A^, A^) of the eigenvalues A = A^ + iA^. However, it is well-known [26] that 
in the standard Hamiltonian case, the relevant “internal mode” of the solitary wave associated with translation has 
a positive energy or signature and hence its collisions with other modes, including ones of the continuous spectrum, 
do not lead to instability. Here, however, as illustrated in Fig. [^exactly the opposite occurs. As the parameter 6/ is 
varied, the relevant eigenvalue moves towards the continuous spectrum (whose lower limit is X = ±i) and the collision 
with it leads to a complex eigenvalue quartet, a feature that would never be possible for a single soliton of the regular 
NFS, under a translation-symmetry-breaking perturbation. This is a remarkable feature of the T^T-symmetric NFS 
model that is worthy of further exploration, possibly utilizing the notion (recently discussed for T^T-symmetric models 
in EH) of Krein signature. Notice that the work of EH considered a case in the vicinity of the T^T-phase transition, 
whereas in our setting, such a phase transition is impossible, given the real nature of the super-partner Poschl-Teller 
potential, as discussed above. 

Fig.|5] shows that dipolar solutions with a single node are unstable for almost all of their existence interval expect 
when 2.43 < bj < 2.5, i.e. in the immediate vicinity of the linear limit. There are three different instability intervals: 
(1) for bj < 0.48 there are two instabilities, one of exponential nature and an oscillatory one, for 0.48 < bj < 1.31 the 
oscillatory instability is the only one that persists, while the formerly real eigenmode crosses the spectral plane origin 
and becomes imaginary for larger bj. (3) for 1.31 < 6/ < 2.43, there are two oscillatory instabilities, the previously 
mentioned one, and another one caused (in a way similar to the nodeless case) by the climbing up the imaginary axis 
of the eigenmode formerly unstable as a real pair, and its eventual collision with an eigenvalue bifurcating from the 
continuous spectrum. 

For nodeless solutions with k = 3, we can observe in Fig. (top panels) that the solution is unstable throughout its 





7 




FIG. 5: Imaginary part (left panels) and real part (right panels) of the eigenvalues associated with the linearization around 
the solution of the nodeless (top panels) and the single-node (i.e., dipole) solutions branches (right panel) for = 1. It can 
be observed that the nodeless solutions become unstable in the interval 0.56 < bj < 1.37, whereas the solutions with a single 
node are unstable for all bi except for a very small interval 2.43 < bi < 2.5 close to the upper existence limit; in addition, for 
bi < 0.48 the latter waveform is also exponentially unstable. 



FIG. 6: Two case examples of the spectral plane (Ar, A^) of eigenvalues A = Ar + iXi of the solution for 5/ = 0 (left panel) and 
5/ = 1 (right panel). The eigenvalue which is spectrally stable in the left panel but whose collision with the band edge of the 
continuous spectrum is responsible for the instability in the right panel is denoted by a red mark. 


range of existence because of an eigenmode entering the phonon band and causing oscillatory instabilities; a second 
localized mode enters at bj = 1.33 and, finally, for 6/ = 1.89, the soliton becomes exponentially unstable. Analogously, 
the solutions with a node are unstable for all 6/, incurring, in fact, typically multiple instabilities for each value of the 
parameter, which can be summarized as follows (see the bottom panels of Fig. an oscillatory instability is present 
for almost every value of bj {bj < 1.82); apart from this we observe, for low values of 6/, two pairs of real eigenvalues 
which coalesce into a quartet at bj ~ 0.525; this quartet leads to two imaginary pairs when bj ~ 1.14; one of them 
moves down along the imaginary axis and finally, at bj ^ 1.27, an additional instability due to a real pair emerges.. 

It is relevant to note in passing, another interesting result which relates to the k < 1 case; we have found that for 
K < 2/3, the nodeless soliton is stable for every bj. This, as well as the results above indicate the strong dependence 
of the stability properties on the precise strength of the nonlinearity parameter. 

Finally, we consider the dynamics of these unstable waveforms for several prototypical cases in Fig.[^ For bj = 0.65 























FIG. 7: Imaginary part (left panels) and real part (right panels) of the eigenvalues associated with the linearization around the 
solution of the nodeless (top panels) and the single-node solution branches (right panel) for = 3. It can be observed that all 
the solutions (for different values of bi) are unstable; see the text for a detailed description of the eigenvalue variation over bj. 


and K = 1 we observe that when t > 200, the oscillatory (as predicted by our eigenvalue computations) nature of 
the instability gradually kicks in and eventually renders the solitary wave more highly localized (i.e., narrower) at 
X = 0 and with a larger amplitude (i.e., taller). Subsequently, the amplitude of the pulse is subject to breathing, 
but it remains fairly robust, even after multiple collisions with small amplitude radiative wavepackets scattering back 
and forth from the boundaries (not visible in the scale of the plot). For 6/ = 1 and k = the growth rate is larger 
and the instability effects stronger; it manifests in an oscillatory growth of the soliton (given the oscillatory nature of 
the instability), as well as a “swinging” of the solution between the gain {x < 0) and loss {x > 0) regions, according 
to the terminology of [24]. This eventually leads to rapid growth (beyond the resolution of the numerical scheme). 
We do not follow the solution past these large values of its amplitude. This behaviour is generic for the oscillatory 
instabilities as long as the growth rate is above a threshold, as shown also in the example for 6/ = 1 and k = 3, and for 
the nodeless and single-node solitons. Finally, we have considered the effect of the exponential instabilities in solitons 
with a node and hz = 1. Those solitons are both exponentially and oscilatorily unstable for bj < 0.48. In that interval, 
the soliton is double-humped (see Fig. [^. In the particular example of Fig. we have taken bj = 0.2 where the 
exponential instability dominates to the oscillatory one. The dynamics here can be described as follows: the hump 
originally located at the loss {x > 0) side shifts to and remains pinned with regular oscillations of the amplitude at 
X = 0. On the other hand, the hump initially at the gain {x < 0) side is “ejected”, as a result of the instability, along 
the (exponentially localized around x = 0) gain side of x < 0. 

IV. CONCLUSIONS & FUTURE CHALLENGES 

In the present work, we revisited a potential that has been explored previously in a number of studies relating to 
T^T-symmetric models. We discussed how for a special monoparametric family within this model, it is not only VT- 
symmetric but also super-symmetric with a partner which is the Poschl-Teller potential, a feature which enabled us to 
identify its purely real spectrum (and the bifurcations of bound states within it) and to corroborate the corresponding 
results numerically. As a byproduct of its super-symmetric origin, this family of potentials was found to be devoid of a 
T^T-pliase-transition. We then turned to a nonlinear variant of the model for arbitrary powers of the nonlinearity, and 
illustrated that exact nonlinear solitonic solutions degenerated in the appropriate limit to the linear states identified 
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FIG. 8: (Top left panel) Space-time contour plot of the evolution of the squared modulus (density) of the solution during its 
unstable dynamics for bi = 0.65 and = 1; the inset shows the evolution of the (maximal) density of the solution occurring 
at X = 0. The top right (bottom left) panel shows the evolution of the density for 5/ = 1 and ac = 1 (/^ = 3); the inset shows 
the evolution of the norm and displays its eventual indefinite growth. The bottom right panel considers the evolution of the 
unstable soliton with a node for bi = 0.2 and = 1, leading eventually to a split of the two humps into a stationary (at x = 0) 
and a traveling one (at x < 0). 


previously. While there was no T^T-phase-transition in this model, we found that the nonlinear solutions were still 
subject to instabilities, such as e.g. the one stemming from a collision of an internal mode with the continuous 
spectrum (band edge), leading to a quartet of eigenvalues. The ensuing oscillatory instability led to an oscillating, 
progressively larger amplitude soliton in the cases examined. Additional families of solutions were discussed, including 
e.g. the one-node branch (dipolar solution), and their reduced stability (in comparison to the nodeless branch) was 
illustrated. 

While this work, to the best of our understanding, is only a first step in connecting all three notions of VT- 
symmetry, super-symmetric potentials and nonlinear phenomenology (including instabilities), naturally this is a 
theme that is worthy of considerable further studies. For one thing, there are numerous additional super-symmetric 
potentials with real spectra that can be devised and are worth examining. For instance, the sl(2, C) considerations 
of [18] already suggest some such options including the super-potentials W{x) = {m — l/2)cothx — i6/cosech(x) 
or W{x) = ±(m — 1/2) — i6/exp(=F^)- Additionally, there have already been proposals for T^T-symmetric square 
well potentials considered in the SUSY framework [28], and for non-Hermitian SUSY hydrogen-like Hamiltonians 
with real spectra [29]. Especially in the latter higher dimensional context, understanding the delicate interplay of 
T^T-symmetry, super-symmetric models with their bound states, and collapse induced by nonlinearity could be an 
especially interesting topic. Finally, there are some potentially intriguing deeper connections. SUSY partners are 
based on commutation formulae as are integrable nonlinear equations. Perhaps the latter is intrinsically responsible 
for the similarity of the structure of the SUSY partner potentials with the well-known Miura transformation respon¬ 
sible for converting the modified Korteweg-de Vries equation to the Korteweg-de Vries equation [30] . Exploring these 
connections further would constitute an important direction for further studies and efforts along this direction are 
already underway m- 
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